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Abstract 

The simulation of held electron emission from arrays of micrometer-long open-ended (5, 5) carbon 
nanotubes is performed in the framework of quantum theory of many electrons. It is found that the 
applied external held is strongly screened when the spacing distance is shorter than the length of 
the carbon nanotubes. The optimal spacing distance is two to three times of the nanotube length, 



slightly depending on the applied external fields. The electric screening can be described by a 
factor that is a exponential function of the ratio of the spacing distance to the length of the carbon 
nanotubes. For a given length, the held enhancement factor decreases sharply as the screening 
factor larger than 0.05. The simulation implies that the thickness of the array should be larger 

o : 

than a value but it does not help the emission much by increasing the thickness a great deal. 

o 

(N 

|>! ■ PACS numbers: 73.22.-f 73.21.-b 79.70.+q 

O. 

C ■ . Keywords: field emission, carbon nanotube, screening effect 

O 

[Xj i I. INTRODUCTION remarkably. 

Q fl s y, fl fl fl B 3 it is 

clear that denser and longer CNT bundle has 
For the application of carbon nanotubes stronger screening effect that would reduce 
(CNTs) as the cold cathodes of field elec- the field electron emission of each CNT. On 
tron emission (FE), the ideal structure is the other hand, larger spacing distance de- 
to arrange the CNTs into an aligned ar- creases the number of CNTs in a unit area, 
ray. It has been observed that the spac- which leads to weaker mean emission cur- 
ing distance of CNTs affects FE properties rent density. It is important to describe the 

screening quantitatively and to find out its 
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effect on the FE ability. The calculation field, thereby it is also related to the screen- 
of Nilsson et al. lj suggested that the opti- ing effect. The enhancement factor of the 
mum spacing distance would be twice the array has been calculated by the classical 
length of the CNTs. However, Jung Sang method. [9j When the charge redistribution 
Suh et a/.[2j showed that the emission cur- and the screening effect are taken into ac- 
rent density is optimized when the length of count, as will be shown in the present paper, 
the CNTs is equal to the spacing distance, the field enhancement factor of the array is 
By solving Laplace's equation, Bocharov and different from the classical one obviously. 
Eletskiijofl found that the emission current To estimate the excess charge distribution 
density has the maximum value when the and the field penetration at the apexes of 
spacing distance is half of the length of the SWCNT arrayS; it requ ires a large-scale sim- 
CNTs. More careful studies on this topic ulation that ghould reflect both the quan _ 
would obviously be useful. In this paper, we tum electron struct ure at the apex and the 
will simulate the FE of the array of single- Coulomb interaction over the tubes. In ex- 
walled carbon nanotubes (SWCNTs) with a perimentS; the length of CNTs is usually in 
quantum/molecular hybrid method. micrometers, while the radius is in nanome- 
Under the applied external field (denoted ters. A huge number of freedoms are in- 
ky F app i), it has been known that the volved. For instance, the (5, 5) type SWCNT 
SWCNT is charged. [1QJ Therefore, the SWC- of 1 jjm length consists of about 10 5 carbon 
NTs of the array are coupled to each other atoms. Limited by the computational effi- 
through the Coulomb interaction. The ex- ciency and resources, all ab initio studies so 
cess charges are also the origin of the screen- far can only simulate the local properties in- 
ing. The competition between the screening volving hundreds of carbon atoms. As the 
effect and the density of CNTs would be com- electronic properties are sensitive to both the 
plicated by the field enhancement factor that detailed atomic arrangement (i.e., the loca- 
ls length-dependent. The field enhancement tion of defects, adsorbates, and the chiral- 
factor was presumably proportional to the as- ity) and the distribution of excess charges 
pect ratio of the SWCNT. However, recent over the whole tube, it is a big challenge 
quantum simulations revealed that the field to simulate a SWCNT array that consists 
penetration at the apex of the SWCNT is of SWCNTs with length in the order of mi- 



significant. 10, The field penetration de- crometer. Only recently has it been possi- 



pends on the value of the applied external ble to tackle an individual SWCNT of re- 



alistic size in the FE conditions by a mul- whose direction is parallel to the axes of the 
tiscale method involving quantum mechan- CNTs, and the S WCNTs are mounted ver- 
ics and molecular mechanics. 
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ll| In the tically and formed a regular square lattice 

present paper, we adopted this method to on the cathode surface. The distance of two 

simulate the FE of the S WCNT array, of nearest neighbored S WCNTs (to be referred 

which the length is in the micrometer scale, to as "spacing distance") and the length of 

The structure of the S WCNTs is specified to SWCNTs in the array are denoted by d and 

the (5, 5) armchair type. The dangling bonds L, respectively. When a field is applied to the 

in the open mouths of the SWCNTs are satu- SWCNT array, the electrons have opportu- 

rated by hydrogen atoms. We further assume nity to emit into vacuum through the apexes 

that the SWCNTs are vertically mounted on of the SWCNTs by quantum tunneling, 

a metal surfaced uniformly. Therefore only g ince electrons are emitted from the apex 

one SWCNT is required to be dealt with in of each SWCNT by quantum tunneling, the 

the environment of other SWCNTs. Each apex part mugt be treated by quantum me _ 

SWCNT of the array is treated as a replica chanicg The part Qn thg gubstrate gide 

of the SWCNT under simulation. mainly affects the field emission through 

In Section II, the simulation method is re- Coulomb potential of the excess charges, so it 



viewed briefly. The simulation results are can be treated by a semiclassical method, [id 

111 ] Therefore we should divide each SWCNT 



presented and discussed in Section III. The 



dependence of the screening on the spacing into a quan t U m region and a semiclassical 

distance, the length, and the applied exter- reg i on . The quantum region is dealt with 

nal fields are discussed in Section IV, where on atom i c sca l e where the density matrix of 

a factor is introduced to describe the screen- tbe e i ec trons is obtained quantum mechani- 

ing, and the correlation between this factor cally The quan t U m region should be large 

and the field enhancement factor is presented. enoug h to ensure that the artificial division 

The last section gives the conclusions. does not affect the physical results seriously. 

By our experience, the proper size of the 

II SIMULATION METHOD quantum region is much bigger than that the 

standard ab initio methods could deal with. 
The CNT array for FE is a typical multi- We have to further divide the quantum re- 
scale system. In our model, the CNT array gion into sub-regions. Each sub-region to- 
is under the uniform applied external field, gether with its adjacent sub-regions forms a 



subsystem that is dealt with by the modi- age charges of the excess charges of the SWC- 
fied neglect of diatomic overlap (MNDO)[l3] NTs. It should be noted that even in the 
semiempirical quantum mechanical method semiclassical region the electron energy band 
(here the MOPAC software has been used), structure originating from the quantum me- 
The excess charges outside the subsystem be- chanics should be taken into account. For 
ing dealt with are treated as point charges, the (5, 5) SWCNT, there are both experi- 
Their contribution to the subsystem being mental and theoretical evidences for the con- 
dealt with is through the Coulomb interac- stant density of state (DOS) in the vicinity 



tion. To accelerate the simulation, we first of the neutrality level. |13l . Il4l . Il5l . Il6| It turns 

simulate an isolated subsystem under vari- out that the excess charge density can be ap- 

ous external fields. Assuming all subsystems proximated by a linear function of the longi- 

except that contains the apex are resemble tudinal coordinate of the SWCNT. [itJ 
to the isolated subsystem, their electron den- The coupling of quantum reg ion and semi- 

sity can be read out from the database, which classical region is through the quasither- 

has been constructed by simulating the iso- modynamic equilibrium condition which as- 

lated subsystem, with the entry of electric gumes that the chemical potentia i ( Fe rmi 

field that is the superposition of the applied level) is a constant over the entire array . To 

external field and the fields contributed by avoid the comp i exity arose from the Schot- 

the excess charges in other parts of the array tky junction that would emerge at the back 

as well as by the image charges. With this ac- contact in principle, we simply assume the 

celeration algorithm, we are able to deal with Fermi leyd of the SWCNTs is 5 eV below 

a quantum region of length over 900 nm. the yacuum potential in the absence of ap . 

Although in principle we can deal with plied external field. The density of the ex- 
an entire isolated SWCNT quantum mechan- cess charge ("excess density" for simplicity) 
ically, the connection of the SWCNT and the calculated separately in quantum and semi- 
cathode is too complicate for a full quantum classical regions should coincide at an over- 
treatment. It is convenient to treat the part lap place of two regions. The self-consistent 
of SWCNT connecting the cathode as the excess density of the entire SWCNT array 
semiclassical region. In the semiclassical re- is achieved through iterations that contain 
gion, the Coulomb potential is governed by a small loop and a big loop. In the small 
Poisson's equation. The boundary condition loop, the subregions of the quantum region 
of the metal surface is guaranteed by the im- are dealt with one by one, and repeated un- 



til that a converged electron density in the E k (ir*)/h, which is approximately equal to 

quantum region is obtained. In the big loop, 10 14 Hz. 

the quantum region and the semiclassical re- We should focus at the electrostatic po- 

gion are dealt with alternatively until the self- tential U(z) in the vicinity of the apex of the 

consistent charge distribution is achieved. SWCNT, which determines the major feature 



of the FE of the array. 



III. SIMULATION RESULTS 



Ti ii i j _.i A. Varying spacing; distance 

It has been shown that the vacuum po- 

tential barrier in the circumambience of the 



SWCNT is high and thick 10], therefore the 
electrons would most probably emit forward 
from the first layer of the tip. The transmis- 
sion coefficient (D) can be estimated by the 
WKB approximation 



D = exp 



'2m[U(z) - E { ]dz 



In Fig. H](a) /(b), we plot U(z) of differ- 
ent spacing distances under the applied ex- 
ternal field, F app i=12.0 V//im. The length of 
the tubes, i.e., the thickness of the array, is 
1.00 /iin in Fig. [T](a) and 0.75 fim in Fig. 
E!_[b). The Z axis has its origin at the last 
atom of the SWCNT and is parallel to the 
direction of the tube axis. The Fermi level 



of the electrons in the nanotubes is assumed 

where Ufz) is the electron energy potential, . , r „ -, r . . , ,. -p, r 

w to be -5.0 ev m our simulation, iherefore 



are m 



Er is the Fermi energy, and the integral is over ,. , n , , . m 

the potential profiles shown m rig. [1] . 

the classical forbidden region where U(z) — r , , , .... ... 

v ' tact the apex-vacuum potential barriers that 



Er > 0. We have assumed that the electrons , n , •, i_i.-Tj.ri i- 

control the probability of electron tunneling. 

possess the Fermi energy. -p, , . c ,, , , . , , , 

Ihe observation of the barrier height iower- 

With D in hand, the emission current of , . 

mg as a increasing is consistent with the b F 

each SWCNT is estimated by , r£ , , , . LI . 

mechanism of field-reduced barrier |10l] since 

larger d should lead to smaller field screen- 

I vq ex cD, (2) increase the effective field applying 

where g exc are the extra electrons of the first to the SWCNTs. The sharper shape of Fig. 

layer atoms, and v is the collision frequency [T^a) comparing with Fig. [T](b) is a conse- 

(the number of electrons hitting the barrier quence of the field enhancement. The curves 

per unit time) that can be estimated from almost coincide with each other for d> 1.5 /.m 

the average kinetic energy of ir* electrons as in Fig. (Ha) /(b), implying that the screening 
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FIG. 1: (Color online) The potential of different 
spacing distances (d) in F app i=12.0 V//im. (a) 
-L=1.00 /im. The curves from up to down are 
corresponding to different d ranged from 0.50 
/im to 2.50 /im with step of 0.25 /im. The po- 
tential of d=4.00 /im and that of the individual 
SWCNT are also plotted, but they are invisible 
since they coincide with the lowest curve, (b) 
L=0.75 /im. The curves from up to down are 
corresponding to different d ranged from 0.50 
/im to 2.50 //m with step of 0.25 /im. The po- 
tential of <i=0.375 /mi and that of the individual 
SWCNT are also plotted, and corresponding to 
the uppermost and lowermost curves. 

effect is negligible when d>1.5 /im. 

The mean current density from the ar- 
ray ("current density" for simplicity and de- 
noted by J) against d are presented in Fig. 
[5] for three sets of parameters. The lengths 
of SWCNTs are 0.75 /im, 1.00 /im, and 1.00 
/im; the external fields are 12.0 V//im, 12.0 
V//im, and 10.0 V//im for the Figs. [2(a), (b), 



and (c), respectively. The current density is 
very sensitive to both L and -F app i (the J axes 
for Figs. E^a), (b), and (c) are multiplied by 
the factors 10~ 5 , 10°, and 10 -3 , respectively). 
It is notable that the emission is turned-on at 
certain spacing distance, which is about 1.0 
/im here, roughly equal to the length of the 
tube. Comparing Figs. [2](a) , (b), and (c), 
one may see that both the turn-on spacing 
distance and the maximum of J would de- 
pend on F app i. 

The parabolic decrease of the current den- 
sity as the spacing distance getting large can 
be easily understood by the fact that the cur- 
rent density is proportional to the number of 
SWCNTs in a unit of area when the screening 
effect is negligible. 

B. Varying length 

We plot the U(z) for different L under the 
applied external field F app i=10.0 V//iin, with 
d=0.75 in Fig. [3](a) and 1.00 /im in Fig. fflb). 
The lowering of apex-vacuum barrier is asso- 
ciated with the excess charge accumulation 
at the apex. For the same L, the barrier of 
larger d is lower than that of shorter d. It 
implies that the array of larger d can accom- 
modate more excess charges at each apex. 

The current density CIS db function of L is 
shown in Fig. |H The squares, circles, and 
triangles in this figure are corresponding to 
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FIG. 2: The current density of the SWCNT 
array versus the spacing distance, (a) L=0.75 
[Ain and F app i=12.0 V//zm. (b) L=1.00 /iin 
and F app i=12.0 V/mn. (c) L=1.00 /mi and 
F appl =10.0 V//im. 




FIG. 3: (Color online) f/(zj for different Z un- 
der the applied external field of 10.0 V//xm. The 
curves from up to down are corresponding to dif- 
ferent L ranged from 0.25 /mi to 4.00 /mi with 
step of 0.25 mm (a) d=0.75 mm (b) (2=1.00 
/mi. 



the data of (2=0.50 /im, (2=0.75 /im, and 
<2=1.00 /im, respectively. From this figure, 
one sees that the emission current densities 
of the array of short SWCNTs are negligi- 
ble. The minimum length for a significant 
current density (say 10~ 10 A /cm 2 ) is enlarged 
as the spacing distance decreases. The rapid 
increase of the current densities is the con- 
sequence of the field enhancement. However, 
when L is large, the screening has significant 
effects and the current densities only increase 
moderately. The array with L=0.5 /im can 
hardly emit electrons. 
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Length (fim) 



FIG. 4: (Color online) The current density ver- 
sus the length of the SWCNTs. The squares, cir- 
cles, and triangles are corresponding to (f=0.50 
jum, d=0.75 /xm, and d=1.00 //m, respectively. 




Z(nm) 

FIG. 5: (Color online) U(z) for the -F a ppl ranged 
from 7 to 15 V/fj,m (up to down) with d=2.00 
/im and L=1.00 /im. 

C. Varying applied external field 

We have considered the array with d=2.00 
//m and L=1.00 /im for various applied ex- 
ternal fields. The electron potential U(z) is 
presented in Fig. [5] for F app i ranged from 7 to 
15 V//xm. 

The J--F app i characteristic is given in Fig. 



FIG. 6: The current density versus the applied 
external field for the S WCNT array with d=2.00 
/im and L=1.00 fim. The insert shows the cor- 
responding Fowler-Nordheim plot. 

|6j The inset of Fig. [6] is the Fowler-Nordheim 
plot. The nonlinear Fowler-Nordheim plot 
implies that the mechanism for the FE of the 
SWCNT array could be different from that 
of the metal plane emitters. 

IV. SCREENING FACTOR 

In order to describe the screening effect 
quantitatively, we define a screening factor a 
as 

V 

a=l- T =^, (3) 

where V is the voltage drop (related to 
the substrate) at the middle point of the 
line connecting two neighborhood apexes. 
It should be zero if there is no screening 
and equal to 1 if the array, as an ideal 
metal layer of thickness L, screens the field 




FIG. 7: (Color online) (a) The screening factor 
a versus the ratio of the spacing distance to the 
length of the SWCNTs (d/L). (b) The screening 
factor a versus F app \. 

completely. The screening factor a versus 
d/L is shown in Fig. [3(a) for four sets 
of parameters: (1) d=0.75 fim, F app i=12.0 
V/fjm] (2) G?=1.00//m, F app i=12.0 V///m; (3) 
L=0.75 /iin, F appl =10.0 V///m; (4) 1=1.00 
jttm, F app i=10.0 V//im. The upper (lower) 
triangles are corresponding to various d with 
fixed L=0.75 (1.00) jum. The squares (circles) 
are corresponding to various L with fixed 
d=0.75 (1.00) fim. Notably, all points fall 
into a curve, implying that the screening fac- 
tor is a function of d/L. The best fitting to 
the simulation (solid curve) is 

a = exp (-(2/0.65L), (4) 

For the array of d=2.00 /im and L=1.00 
jttm, the screening factor a is calculated with 
different F appl . As shown in Fig. [7(b), we find 
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that it is almost independent of the applied 
external field. It implies that the screening 
factor fixed by the ratio d/L is an intrinsic 
feature of the array. 

From the data of U(z), we can estimated 
the local field strength (F m ) as the value of 
the slope of the barrier potential U(z) at the 
point of steepest decrease. An effective field 
enhancement factor f3 can then be defined as 
F m /F app i. Figure [8] shows the correlation be- 
tween f3 and the tube length (a), and the cor- 
relation between (3 and the screening factor 
(b), for fixed F app i and d . The field enhance- 
ment factor of the array is not proportional 
to the tube length, which is completely differ- 
ent from an individual nanotube (or a metal 
tip). The increase of tube length strengthens 
(3 and the screening tends to reduce (3 . But 
the increase L will lead to increase of screen- 
ing. Therefore, the dependence of j3 on L and 
a is complicate. If there were no screening, 
(3 in Fig. [5(b) should increase faster. 

The reduction of field enhancement by 
screening can be seen more clearly in Fig. 
for a few sets of applied external fields and 
tube lengths. Figure [9(a) shows that for large 
spacing distance, where the screening is week, 
the field enhancement factor tends to a con- 
stant. For a given L, the classical field en- 
hancement factor should be a constant. How- 
ever, as one seen in [9(b), that would be true 
only if the screening could be ignored. For 




Length (|im) a 

(a) (b) 

FIG. 8: (Color online) The field enhancement 
factor (/?) in F app i=10.0 V//iin. The squares are 
for d=0.75 /im, and the triangles are for d=1.00 
/iin. (a) (3 versus the length of the SWCNTs; (b) 
(5 versus the screening factor. 

the instances of the [2(b), the field enhance- 
ment factor decreases sharply as a increase 
when 1/a < 20. 

V. CONCLUSIONS 

We have studied the field electron emis- 
sion properties of arrays of single-walled 
carbon nanotubes in the quantum level. 
As the spacing distance of the SWCNTs 
exceeds a certain value (about the length 
of the tube), the emission current density 
increases rapidly. When the spacing distance 
is large enough, the current density decreases 
with the spacing distance parabolic 
function. An optimal spacing distance that 
corresponds to the maximum current density 
could be two to three times of the nanotube 
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Distance (urn) 1/a 

(a) (b) 

FIG. 9: (Color online) (a) The field enhance- 
ment factor (/?) versus the spacing distance of the 
SWCNTs (d); (b) (5 versus 1/a, for three sets of 
parameters: (1) squares: -F app i=12.0 V//im and 
L=0.75 /im, (2) triangles: F app i=10.0 V//im and 
L=1.00 /im, (3) circles: F app i=12.0 V//«n and 
L=1.00 /im. 

length. The exact value of the ratio corre- 
sponding to the peak current density would 
depend on the applied external field and the 
tube structure. That would be one of the 
reasons for the discrepancy of experimental 
results of different groups. The screening 
effect can be described by the factor (a) 
denned in Eq. [31 which is found to be an 
exponential decreasing function of the ratio 
of the spacing distance to the length of 
the SWCNTs, and is almost independent 
of the applied external field. Therefore, 
the screening factor introduced here reflects 
an intrinsic character of the array. When 
a < 0.05, the screening effect can be ignored 



and each SWCNT of the array behaves as an 
individual emitter. For large a, on the other 
hand, the emission property of a SWCNT 
in the array is obviously different from the 
individual SWCNT. The field enhancement 
factor of the SWCNT array is not a linear 
function of the tube lengths as the metal 
rod model predicted. And it depends on 
the applied external fields and the spacing 
distances. For a given tube length, it is a 
monotonically decreasing function of the 
screening factor. Although the field enhance- 
ment factors depends on the spacing distance 
(Fig. H](a)) in the same trend as the classical 
calculation [9], the quantitative discrepancy 
between quantum simulation and classical 
calculation is large. To increase the emission 
current density, one can either adjust the 
spacing distance or increase the lengths of 
the tubes. For the (5,5) SWCNT array 
considered in the present paper, we find 
that the current density is very small for the 
spacing distance 0.5 /xm. For larger spacing 
distance (0.75 and 1.0/xm, for instances), 
the current density increases rapidly as the 
length increases until the length is as large as 
1.5 times of the spacing distance. For longer 
length, the current density only increases 
slowly. It implies that lengths of SWCNTs 
of the arrays need not be too long. 
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